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We report searches for B-meson decays to the charmless final states pK* and fo(980)K* with 
a sample of 232 million BB pairs collected with the BABAR, detector at the PEP-II asymmetric- 
energy e + e~ collider at SLAC. We measure the following branching fractions in units of 10 -6 : 
B{B+ -> p°K *+) = 3.6 ± 1.7 ± 0.8 (< 6.1), B(B + -> p + K*°) = 9.6 ± 1.7 ± 1.5, B(B° -> p~K* + ) = 
5.4 ± 3.6 ± 1.6 (< 12.0), B(B° -> p° K*°) = 5.6 ± 0.9 ± 1.3, B(B + -» / (980)iT + ) = 5.2 ± 1.2 ± 0.5, 
and B(B° -> f (980)K*°) = 2.6 ± 0.6 ± 0.9 (< 4.3). The first error quoted is statistical, the 
second systematic, and the upper limits, in parentheses, are given at the 90% confidence level. For 
the statistically significant modes we also measure the fraction of longitudinal polarization and the 
charge asymmetry: f L {B + -» p + K*°) = 0.52 ± 0.10 ± 0.04, f L (B° -> p°K*°) = 0.57 ± 0.09 ± 0.08, 
Acp(B + -» p+K*°) = -0.01 ± 0.16 ± 0.02, Acv(B° -> p°K*°) = 0.09 ± 0.19 ± 0.02, Acp(B + -> 
/o(980)ff* + ) = -0.34 ± 0.21 ± 0.03, and A C p(B a -> f (98Q)K*°) = -0.17 ± 0.28 ± 0.02. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



The study of B-meson decays to charmless hadronic 
final states plays an important role in understanding 
CP violation. The charmless decays B —* pK* pro- 
ceed through dominant penguin loops and Cabibbo- 
suppressed tree processes (B+ — ► p + K*° is pure pen- 
guin) to two vector particles (VV). A large longitudinal 
polarization fraction fa (of order (1 — 4m 2 / /m 2 3 ) ~ 0.9) 
is predicted for both tree and penguin dominated VV 
decays p|. However, recent measurements of the pure 
penguin VV decays B — > <f>K* indicate fa ~ 0.5 0- 
Several attempts to understand this small value of fa 
within or beyond the Standard Model (SM) have been 
made Q. Further information about S'f7(3)-related de- 
cays may provide some insight into this polarization puz- 
zle. Characterization of the four B — > pK* modes can 
also be used within the SM framework to help constrain 
the angles a and 7 of the Unitarity Triangle [4j . 

We report measurements of branching fractions, longi- 
tudinal polarizations, and direct CP-violating asymme- 
tries for the B — » pK* decay modes. We also measure 
branching fractions and direct CP- violating asymmetries 
for the B — > /o(980)A"* modes that share the same fi- 
nal states. We present improved analyses of previously 
measured modes with larger statistics and explicit 
consideration of non- resonant backgrounds. We measure 
B° -► p-K* + , B° -> p°K*°, B+ -> f (980)K* + , and 
B° -> /o(980)ir" for the first time. Charge-conjugate 
modes are implied throughout. 

This analysis is based on a data sample of 232 million 
BB pairs, corresponding to an integrated luminosity of 
210 fb" 1 , collected with the BABAR detector @ at the 
SLAC PEP-II asymmetric-energy e + e~ collider operat- 
ing at a center-of-mass (cm.) energy y/s = 10.58 GeV, 
corresponding to the T(45) resonance mass. 

The angular distribution of the pK* decay products, 
after integrating over the angle between the decay planes 
of the vector mesons, for which the acceptance is uniform, 
is proportional to 

i(l - fa) sin 2 8 K > sin 2 6 p + fa cos 2 8 K * cos 2 8 p , (1) 



where 9k * and 9 p are the helicity angles of the K* and 
p, defined between the K(tt + ) momentum and the di- 
rection opposite to the B in the K*(p) rest frame Q- 
We also measure the time-integrated direct CP- violating 
asymmetry _4 C p = (r~ - r+)/(r _ + T+), where the su- 
perscript on the total width T indicates the sign of the 
6-quark charge in the B meson. 

We fully reconstruct charged and neutral decay prod- 
ucts including the intermediate states p° or /o(980) — > 

7T+7T-, p+ -> 7T+7T , K*° -> K+TT~ , K* + -> K+TT°, 

K* + -> K^tt+ (only in p°K* + ), tt° -> 77, and K° s -> 
7T+7T - . We assume the /o(980) measured lineshape || 
and a branching ratio of 100% for /q(980) — > 7r + 7r _ . Ta- 
blen]lists the selection requirements on the invariant mass 
and helicity angle of B-daughter resonances. 

TABLE I: Selection requirements on the invariant mass (in 
GeV) and helicity angle of B-daughter resonances. 

Mode m-n-n mn^ cos 6 P cos 9 k * 

p u A"* + + 7r0 (0.52,1.10) (0.75,1.05) (-0.95,0.95) (-0.5,1.0) 
p K*+ l + (0.52,1.10) (0.75,1.05) (-0.95,0.95) (-0.9,1.0) 

p+K*° (0.40,1.15) (0.77,1.02) (-0.66,0.95) (-0.95,1.0) 
P~K*+ +ir0 (0.40,1.15) (0.77,1.02) (-0.80,0.98) (-0.80,0.98) 
p o K *<r * (p.52,1.15) (Q.77,1.02) (-0.95,0.95) (-0.95,1.0) 



Charged tracks from the B-meson candidate are re- 
quired to originate from the interaction point. Looser 
criteria are applied to tracks forming Kg candidates, for 
which we require \m^+^- — rn K o\ < 12MeV, a mea- 
sured proper decay time greater than five times its un- 
certainty, and the cosine of the angle between the re- 
constructed flight and momentum directions to exceed 
0.995. Charged particle identification provides discrimi- 
nation between kaons and pions, and is also used to reject 
electrons and protons. We reconstruct ir° mesons from 
pairs of photons, each with a minimum energy of 30 MeV 
(p°K* + ) or 50 MeV (p+K*° and p~K* + ). The invariant 
mass of 7T° candidates is required to be within 15 MeV 
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(p°K* + ) or 25MeV {p+K*° and p~K* + ) of the nominal 
mass 9] . 

B-meson candidates are characterized kinematically by 
the energy difference AE = E* B — y/s/2 and the energy- 

1/2 

substituted mass mss = [( s /2 + • \>b) 2 /Ef — p^] , 
where (-Bj,pj) and (Eb,Pb) are the four-momenta of 
the T(45) and fl candidate respectively, and the aster- 
isk denotes the T(4S) frame. Our signal lies in the re- 
gion \AE\ < 0.1 GeV and 5.27 < m ES < 5.29 GeV. 
Sidebands in toes and AE are used to characterize the 
continuum background. The average number of signal 
B candidates per selected data event ranges from 1.05 
to 1.27, depending on the final state. A single candi- 
date per event is chosen as the one with the smallest 
B vertex-fit x 2 (p + if*°and p°K*°), the smallest value 
of x 2 constructed from deviations of reconstructed n° 
masses from the expected value (p~K* + ), or randomly 
(p°K* + ). Monte Carlo (MC) simulation shows that up to 
38% (23%) of longitudinally (transversely) polarized sig- 
nal events are misreconstructed with one or more tracks 
originating from the other B in the event. 

To reject the dominant qq continuum background we 
require |cos0t| < 0.8, where 8t is the cm. frame an- 
gle between the thrust axes of the -B-candidate and that 
formed from the other tracks and neutral clusters in 
the event. We also use as discriminant variables the 
polar angles of the B-momentum vector and the B- 
candidate thrust axis with respect to the beam axis, 
and the two Legendre moments Lq and Li of the en- 
ergy flow around the i?-candidate thrust axis in the cm. 
frame jlfll | . These variables are combined in a Fisher dis- 
criminant T (p°K* + ) or a neural network (NN) (other 
modes). Finally, we suppress background from B decays 
to charmed states by removing signal candidates that 
have decay products consistent with D° — > K~ir + (ir°) 
and D~ — > K + tt~i:~ decays. 

We use an extended (not extended in the p + K *°mode) 
unbinned maximum-likelihood (ML) fit to extract signal 
yields, asymmetries, and angular polarizations simulta- 
neously. We define the likelihood for each event candi- 
date i as the sum of rijVj (xi ; a) over hypotheses j (signal, 
qq background, and several BB backgrounds discussed 
below), where the Vj(xi',a) are the probability density 
functions (PDFs) for the measured variables Xi, and rij 
are the yields for the different hypotheses. The quanti- 
ties a represent parameters in the expected distributions 
of the measured variables for each hypothesis. They are 
extracted from MC simulation and (m,ES,AE) sideband 
data. They are fixed in the fit except for some shape 
parameters of the continuum AE and m^s distributions. 
The extended likelihood function for a sample of N can- 



didates is L = exp (— Y^, n j) Ili=i £*• 

The fit input variables Xi are i«es, AE, NN or J 7 , 
invariant masses of the candidates p (/o(980)) and K* , 
and helicity angles 9 P and 9k* ■ We study large control 



samples of B — > Dtt decays of similar topology to verify 
the simulated resolutions in AE and tti-es, adjusting the 
PDFs to account for any difference found. 

Since almost all correlations among the fit input vari- 
ables are found to be small, we take each Vj to be the 
product of the PDFs for the separate variables with the 
following exceptions where we explicitly account for cor- 
relations: the correlation between the two helicity angles 
in signal, the correlation due to misreconstructed events 
in signal, and the correlation between mass and helic- 
ity in backgrounds. The effect of neglecting other cor- 
relations is evaluated by fitting ensembles of simulated 
experiments in which we embed the expected numbers 
of signal and charmless B-background events, randomly 
extracted from fully-simulated MC samples. 

We use MC-simulated events to study backgrounds 
from other B decays. Charmless £>-backgrounds are 
grouped into up to 11 classes with similar topologies 
depending on the mode. Yields for decays with poorly 
known branching fractions are varied in the fit with those 
remaining kept fixed to their measured values. One to 
four additional classes account for neutral and charged B 
decays to final states with charm. Up to 6 classes account 
for misreconstructed events in signal. We also introduce 
components for non-resonant backgrounds such as ixixK* , 
pKir, f (980)Kir, and / (1370)#7T, which differ from sig- 
nal only in resonance mass and helicity distributions. 
The magnitudes of these components are determined by 
extrapolating from fits performed on a wider mass range 
reaching to higher mass values and are fixed in the fit. 
Fig- El shows the sPlots for the invariant mass of Kit 
and 7T7T in the p + K *°and p°K *°modes, respectively. The 
data events are weighted by their probability to be sig- 
nal, calculated from the signal and backgrounds PDFs of 
the AE, mES, and NN variables. 
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FIG. 1: sPlots for the invariant mass of Kir in p + K*° (left) 
and 7T7T in p°if*°and f (980)K*° (right) up to the higher- 
mass regions. The points with error bars show the data, and 
the solid (dashed) lines show the projected PDFs of the sig- 
nal and non-resonant background (non-resonant background 
only: pKn in p + K*°; the sum of fo(1370)K* , ttttK*, and 
tvtyKtt in p°K*°). The arrows show the nominal fit regions. 

The results of the ML fits are summarized in Table ITTI 
For the branching fractions, we assume equal production 
rates of B + B~ and B°B°. The significance S of a sig- 
nal is defined by Aln£ = S 2 /2, where Aln£ represents 
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the change in likelihood from the maximal value when 
the number of signal events is set to zero, corrected for 
the systematic error defined below. We find significant 
signals for p+K*°, p°K* a , and f Q {980)K* + , and some ev- 
idence for /o(980)isr* . For the modes with significance 
smaller than five standard deviations we also measure 
the 90% confidence level (C.L.) upper limit, taking into 
account the systematic uncertainty. Fig.[3shows projec- 
tions of the fits onto tties- 




5.25 5.26 5.27 5.28 5.26 5.27 5.28 5.29 

m ES (GeV) 

FIG. 2: Projections of the multidimensional fit onto rriES for 
events passing a signal-to-total likelihood probability ratio cut 
with the plotted variable excluded for (a) p°K* + , (b) p + K*° , 
(c) p-K* + , (d) p°K*°, (e) / (980)iT + , and (f) / (980)iT . 
The points with error bars show the data; the solid, dashed 
and dotted lines show the total, background, and continuum 
PDF projections respectively. 

A source of systematic error is related to the determi- 
nation of the PDFs and is due to the limited statistics 
of the Monte-Carlo and to the uncertainty on the PDF 
shapes. We obtain variations in the yields ranging from 
1 to 18%, depending on the mode. The systematic error 
due to the non-resonant background extrapolation and 
interference with signal is in the range 6-21%. Event 
yields for B-background modes fixed in the fit are varied 
by their respective uncertainties. This results in a sys- 
tematic uncertainty of 2-12%. We evaluate and correct 
for possible fit biases with MC experiments. We assign a 
systematic uncertainty of 1-7% for this. 

The reconstruction efficiency depends on the decay po- 
larization. For the p°K* + mode we calculate the effi- 
ciency using the measured polarization (combined for the 
two p°K* + modes) and assign a systematic uncertainty 
corresponding to the total polarization measurement er- 



ror (9 and 20% for each mode respectively). For the other 
modes we exploit the correlation between B and fi, and 
obtain the values of B from fits where B and /l are free 
parameters. Fig. |31 shows the behavior of — 21ii£(/l,S) 
for the modes with significant signal. 




BR(B 



V) den 



K*°p°) (10*) 



FIG. 3: Distribution of -21n£(B,/ L ) for B+ -> p+K*°{leit) 
and B° — > *°(right) decays. The solid dots corre- 
spo nd to the centr al values and the curves give contours in 
A x /-21n£(B,/ L ) = 1 steps. 

Additional reconstruction efficiency uncertainties arise 
from tracking (3-5%), particle identification (1-2%), ver- 
tex probability (2%), track multiplicity (1%) and thrust 
angle (1%). Kg and tt° reconstruction contribute 2.3% 
and 3% uncertainty, respectively. Other minor system- 
atic effects are from uncertainty in daughter branch- 
ing fractions, MC samples statistics, and number of B 
mesons. The absolute systematic uncertainty in /x, takes 
into account PDF shape variations (5-10%), B and non- 
resonant backgrounds (4-8%), and efficiency dependence 
on the polarization (1-2%). The absolute uncertainty 
in the charge asymmetry due to track charge bias is less 
than 1%. PDF variations and fixed £>-background effects 
contribute up to 2%. 

In summary, we have searched for B — > pK* and B — > 
/ (980)K* decays. We observe B+ -> p+K*°, B° -> 
p"K*°, B+ -» / (980)iY*+, and B° -> f (980)K*° with 
7.1, 5.3, 5.0, and 3.5 a significance respectively. We mea- 
sure the branching fractions or 90% C.L. upper limits, 
the fractions of longitudinal polarization, and the charge 
asymmetries, summarized in Table ITT1 The measured po- 
larization in the p + K*° and p°K*° modes agrees with 
values measured in (j>K* decays. 
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TABLE II: Summary of results for the measured B-decay modes: signal yield n S i g and its statistical uncertainty, reconstruction 
efficiency e, daughter branching fraction product significance S (systematic uncertainties included), measured branching 

fraction B, (90% C.L. upper limit in parentheses), measured longitudinal polarization }l (for the modes with non-significant 
signals the numbers, in brackets, are not quoted as measurements) and charge asymmetry .Acp- 



Mode 




e(%) Il^(%) S(a) 


S(KT b ) 


h 


-4cp 


p"K* + 




7.9 


32.9 


2.5 
1.3 


3.6 + l;6 ±0.8 (6.1) 
3.2±|.I ± 0.9 


[0.9 ± 0.2] 

[o.s±C:I] 






32t\ 9 r 


15.8 


22.9 


2.1 


3.8^1 ±0.9 


[1.0 ±0.3] 




p+K*" 


f94 ± 29 


13.5 


66.7 


7.1 


9.6 ±1.7 ±1.5 


0.52 ±0.10 ±0.04 


-0.01 ±0.16 ±0.02 




60+£ 


15.2 


32.5 


1.6 


5.4+^ ± 1.6 (12.0) 


[-o.is + r 7 4] 




p«K J 


f85 ± 30 


22.9 


66.7 


5.3 


5.6 ±0.9 ±1.3 


0.57 ±0.09 ±0.08 


0.09 ±0.19 ±0.02 


/o(980)/r + 








5.0 


5.2 ±1.2 ±0.5 




-0.34 ±0.21 ±0.03 


^/o(980)^* + + ^ 




8.5 


32.9 


3.8 


6.2+2-1 ±0.7 




-0.50 ±0.29 ±0.03 


^/o(980)/T+o T+ 


37±n 


16.6 


22.9 


3.2 


4.2±1;| ± 0.5 




-0.13 ±0.30 ±0.01 


/o(980)/-T u 


83 ± 19 


21.7 


66.7 


3.5 


2.6 ±0.6 ±0.9 (4.3) 




-0.17 ±0.28 ±0.02 
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Measurements of branching fractions, polarizations, and direct CP-violation 
asymmetries in B — » pK* and B — > f (980)K* decays 

We report searches for B-meson decays to the charmless final states pK* and /o(980)if* with 
a sample of 232 million BB pairs collected with the BABAR, detector at the PEP-II asymmetric- 
energy e + e~~ collider at SLAC. We measure the following branching fractions in units of 10~ 6 : 
B(B+ -> p°K*+) = 3.6 ± 1.7 ± 0.8 (< 6.1), B{B+ -> p+K*°) = 9.6 ± 1.7 ± 1.5, B(B° -> p~i^* + ) = 
5.4 ± 3.6 ± 1.6 (< 12.0), B(B° -> p°A"* ) = 5.6 ± 0.9 ± 1.3, B(B + -> / (980)ir + ) = 5.2 ± 1.2 ± 0.5, 
and B{B° -> / (980)A"*°) = 2.6 ± 0.6 ± 0.9 (< 4.3). The first error quoted is statistical, the 
second systematic, and the upper limits, in parentheses, are given at the 90% confidence level. For 
the statistically significant modes we also measure the fraction of longitudinal polarization and the 
charge asymmetry: f L {B + -> p + K*°) = 0.52 ± 0.10 ± 0.04, f L (B° -> p°K*°) = 0.57 ± 0.09 ± 0.08, 
Acp(B + -> p + K*°) = -0.01 ± 0.16 ± 0.02, Acp(B° -> p°K*°) = 0.09 ± 0.19 ± 0.02, Acp(B + -» 
/o(980)A"* + ) = -0.34 ± 0.21 ± 0.03, and .4 C p(B° -> /o(980)^*°) = -0.17 ± 0.28 ± 0.02. 



